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摘要 :针对 我 国 海洋 二 号 卫星 微波 散射 计 (HaiYang-2 Scatterometer, HY2-SCAT) MH GEN Æ 
间 分 辨 率 较 低 .无 法 满足 对 于 地 球 冰 冻 圈 .生物 圈 探 测 需 求 的 问题 ,以 其 数据 拓展 应 用 为 目标 ,开展 
了 微波 散射 计 高 分 辨 率 后 向 散射 系数 (co" ) 图 像 重 构 技 术 仿 真 研究 。 采 用 仿真 的 方法 实现 并 对 比 了 
AART(Additive Algebraic Reconstruction Technique), MART (Multiplicative Algebraic Recon- 
struction Technique) ,SIR(Scatterometer Image Reconstruction) 3 种 分 辨 率 增强 的 图 像 重 构 算法 ， 
得 出 了 AART 算法 没有 噪声 抑制 能 力 `\MART 算法 噪声 抑制 能 力 较 弱 、SIR 算法 具有 较 强 的 噪声 
抑制 能 力 , 能 更 好 地 重 构图 像 ,达到 分 辩 率 增强 效果 的 结论 。 然 后 采用 HY2-SCAT 测量 LIB 数据 
对 结论 进行 了 验证 , 重 构 了 岛屿 的 地 貌 特 征 , 得 到 了 分 辩 率 增强 、 细 节 丰 富 、 边 界 清晰 的 重 构 图 像 ; 
重 构 后 的 图 像 还 可 以 应 用 于 植被 . 冰 盖 和 台风 观测 等 地 球 地 理学 研究 领域 。 

关 E 词 :海洋 二 号 卫星 ;微波 散射 计 ; 图 像 重 构 ;分 辨 率 增强 

中 图 分 类 号 :TP 732.1 ”文献 标志 码 :A 文章 编号 :1004-0323(2015)03-0495-10 


术 等 ,但 是 这 些 方法 均 受 限 于 天 线 的 尺寸 。David 
So) fet TAA AART 算法 .MART 算法 以 及 
SIR 算法 对 散射 计 测 量 数 据 进 行 图 像 重 构 ,增强 微 
波 散 射 计 的 测量 数据 分 辨 率 。 这 3 种 算法 在 散射 计 
图 像 重 构 和 分 辩 率 提高 方面 做 出 了 重要 贡献 ,被 广 
泛 应 用 于 陆地 探测 .植被 研究 和 极地 海 冰 探测 。 
海洋 二 号 卫星 散射 计 (HaiYang-2 Scatterome 
ter, HY2-SCAT) 是 我 国 第 一 个 可 业务 化 运行 的 星 
载 微波 散射 计 。 它 的 固有 空间 分 辨 率 大 约 为 30 一 
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随 着 海洋 观测 重要 性 的 进一步 提高 , 星 载 微波 
散射 计 成 为 海洋 动力 环境 观测 的 重要 遥感 器 " ” 。 星 
载 微波 散射 计 因 具有 观测 刘 幅 宽 ,能 在 短 时 间 内 对 
全 球 进行 覆盖 并 提供 高 精度 的 后 向 散射 系数 (o" ) 观 
测 数据 等 特点 ,被 广泛 应 用 于 海洋 风 场 海浪 场 、 数 
值 天 气 预报 以 及 海洋 和 大 气 研究 ”。 由 于 传统 星 载 
微波 散射 计 分 辩 率 较 低 ,通常 在 几 十 干 米 量 级 ,无 法 
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满足 一 些 如 近 岸 风 场 等 千 米 量 级 的 海洋 现象 的 测 
量 ; 同 时 , 随 着 陆地 探测 需求 的 增加 ,迫切 需要 星 载 
微波 散射 计 应 用 于 植被 . 积 雪 、 冰 川 和 土壤 等 的 监测 
和 探测 研究 。 由 于 上 述 拓展 应 用 需求 ,提高 散射 计 
分 辨 率 具有 非常 重要 的 意义 。 通 常 提高 散射 计 分 辩 
率 的 方法 有 多 普 勒 滤波 、 距 离 门 滤波 和 脉冲 压缩 技 
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40 km ,提供 空间 网 格 25 kmX25 km 的 后 向 散射 数 
据 产品 用 于 海面 风 矢 量 场 的 反 演 。 该 散射 数据 产品 
用 于 大 陆 和 海 冰 的 观测 应 用 时 ,分 辨 率 较 低 的 缺点 
更 为 突出 。 为 了 使 HY2-SCAT 能 够 进行 更 广泛 的 
应 用 ,需要 对 其 数据 进行 图 像 重 构 ,获得 更 高 的 分 辩 
率 。 本 文 进行 了 7 AART、MART.SIR 3 种 算法 用 于 
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星 载 微 波 散射 计 图 像 分 辨 率 增强 算法 的 研究 ,并 用 
仿真 数据 和 HY2-SCAT 测量 数据 进行 验证 。 


2 星 载 微波 散射 计 图 像 重 构 模 型 与 方法 


星 载 微波 散射 计 对 地 面目 标 进行 观测 时 ,其 测量 
目标 的 分 辩 率 由 散射 计 天 线 孔径 决定 ,得 到 的 后 向 散 
HRAC ) 是 测量 目标 单元 内 测量 值 的 平均 值 ,导致 
测量 值 和 真实 信号 之 间 存 在 误差 ,因此 需要 用 图 像 重 
构 方 法 对 目标 进行 重 构 。 散 射 计 对 同一 区 域 的 多 次 、 
独立 的 测量 之 间 存 在 重 晋 区域 ,可 以 利用 这 些 重 二 区 
域 的 信息 进行 散射 计 测量 目标 图 像 重 构 "。 

假设 测量 目标 图 像 f 在 有 限 尺 寸 上 足够 连续 ， 
可 以 认为 该 图 像 在 这 个 尺度 上 是 有 限 采 样 或 者 离散 
的 .因此 可 以 把 图 像 f 认为 是 由 比 测量 采样 小 的 、 离 
散 的 、 大 小 一 致 的 像 元 组 成 的 ,所 以 ,第 j 次 无 噪声 
的 测量 值 > 包含 了 很 多 小 的 像 元 , 即 : 

1 M 
z; = g aa has: (1) 
其 中 :* 表示 图 像 了 的 像 元 按 行 扫描 得 出 的 向 量 ; 的 
元 素 , 表示 第 j 次 测量 、 在 第 i 个 像 元 上 的 有 效 孔 
FRAAIE. M 表示 图 像 中 像 元 的 个 数 ,NN 表示 测 
SAAS. Bh OREN TR. 
2.1 AART 算法 

AART 迭代 算法 首先 对 图 像 中 的 每 个 像 元 赋 初 
值 S ( 初 值 通常 是 一 个 相同 的 数 ) ,然后 通过 加 性 更 新 
项 进行 运 代 更 新 ,使 其 不 断 逼 近 真 实 值 .每 个 像 元 的 
更 新 项 是 测量 值 和 前 向 估计 值 的 差 值 ,然后 将 更 新 项 
与 当前 估计 值 求 和 ,从 而 实现 对 像 元 的 更 新 “…“]。 

stl = gf +5 Doha ey — ft) (2) 
其 中 :M 表示 图 像 中 像 元 的 个 数 ; N 表示 测量 的 次 
BD h PERERA TS Ms, ABR 次 迭代 对 > 的 估计 
值 ,p; 表示 覆盖 到 像 元 i 的 测量 的 总 次 数 , fi 为 前 向 
上 映射, fa Fp, 分别 为 : 


I% 
f= — > hs (3) 
qi 1=1 
N 
pi = Dh (4) 
l=1 
其 中 :9 表示 第 7 次 测量 覆盖 的 像 元 的 个 数 , 即 为 : 
M 
qq; = Shi (5) 
i=1 


2.2 MART 算法 
MART 算 法 与 AART 算 法 类 似 , 通 过 迭代 算法 
对 图 中 的 每 一 个 像素 点 进行 更 新 ,与 AART 不 同 的 
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是 , 它 的 迭代 更 新 项 是 乘 性 的 ,每 一 个 像 元 的 更 新 项 
是 测量 值 与 前 向 估计 值 的 比值 ,再 与 当前 估计 值 相 
乘 , 完 成 对 像 元 的 更 新 。 

EB k TORRE. RS 的 每 一 个 像 元 sx 
用 下 列 更 新 项 进行 更 新 : 

sith = sf oe (F) (6) 

其 中 : 参数 w 为 可 调 人 参数 ,w = 1 是 无 加 权 的 
MART,w 关 1 是 加 权 的 MART。 应 用 于 有 噪声 的 
散射 计数 据 时 ,根据 经 验 取 w = 1/2 时 可 以 得 到 较 
好 的 重 构 图 像 ,为 了 简便 迭代 式 写成 如 下 形式 : 


N 
de (7) 
Pi j= 
其 中 ,更 新 项 wj 如 下 定义 : 
us = (a = sid} (8) 
Ho, 
dt = (F) (9) 
2.3 SIR 算法 


当 测 量 图 像 的 信 噪 较 低 时 ,MART 算法 因 无 法 
最 终 收 敛 ,从 而 无 法 取得 较 好 的 图 像 重 构 效 果 . 针 对 
散射 计 目标 测量 的 信 噪 比较 低 的 特点 ,David? fe 
出 了 对 MART 迭代 算法 改进 的 SIR 算法 ,增强 了 算 
法 的 噪声 容 限 ,得 到 了 更 好 的 重 构 效 果 . 改 进 的 算法 
用 一 种 非 线 性 的 权重 和 前 向 测量 估计 代替 了 
MART 算法 中 的 wi ,具体 公式 如 下 : 

p 

[oraa o>) 


[7402 +sta] dt <1 


SIR 算法 通过 对 比例 因子 dj ( 即 测量 值 与 前 向 
估计 的 比值 的 平方 根 ) 大 小 的 判断 ,而 改变 更 新 项 。 
在 理想 环境 下 , 随 着 迭代 次 数 的 增加 ,比例 因子 会 达 
到 单位 值 1, 说 明 所 有 像 元 的 前 向 估计 和 测量 值 相 
一 致 ,迭代 收敛 。 但 是 在 散射 计 实 际 测 量 情况 下 , 因 
为 存在 噪声 ,比例 因子 不 会 最 终 收敛 到 1, 由 于 SIR 
算法 根据 比例 因子 的 大 小 提供 了 不 同 的 非 线性 更 新 
项 ,这 就 限制 了 对 前 向 估计 和 测量 值 之 间 的 大 的 偏 
差 的 更 新 ,从 而 限制 了 最 终 估计 值 的 大 小 ,对 噪声 起 
到 了 限制 作用 ,降低 了 算法 对 噪声 的 敏感 程度 。 


3 重 构 算法 仿真 


3.1 HY2-SCAT 系统 参数 与 性 能 指标 
海洋 二 号 卫星 是 我 国 第 一 颗 海 洋 动力 环境 卫 
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星 ,主要 用 于 全 球 海面 风 场 观测 。HY2-SCAT 采用 
笔 形 波 束 圆 锥 扫描 体制 ,采用 VV( 外 波束 ) .HH( 内 
波束 ) 两 种 极 化 方式 ,通过 天 线 圆锥 扫描 和 卫星 运 
动 , 对 同一 目标 区 域 进行 外 波束 前 视 、 内 波束 前 视 、 
内 波束 后 视 和 外 波束 后 视 4 次 扫描 ,分 别 得 到 测量 
值 a。 海洋 二 号 散射 计 的 内 波束 入 射 角 为 41 ,外 
波束 入 射 角 为 48" ;内 波束 3 dB 波束 宽度 为 
1.5"X1.5°, 外 波束 3 dB 波束 宽度 为 1. 5 X1.5 ;所 对 
应 的 内 波束 足 印 大 小 为 34 kmX26 km, 外 波束 足 印 
大 小 为 42 kmX28 km”, HBS RIN AAR 
样 间隔 决定 ,方位 向 上 的 采样 间隔 由 卫星 地 面 速度 和 
天 线 旋 转速 度 来 决定 ,距离 向 的 采样 间隔 由 脉冲 重复 
频率 (PRF) 和 天 线 旋转 速度 决定 。 卫 星 地 面 速度 约 
为 6. 4 kmy/s, 天 线 旋转 速度 为 16. 67 r/min, 距 离 向 
地 面 足迹 间隔 为 23 km, 脉 冲 重复 频率 为 185 Hz, 
方位 向 内 波束 地 面 足迹 间隔 为 13 km, 外 波束 足迹 
间隔 为 18 km, BSRMBDA 1/3 PKA” , 保 
证 了 重 构 算法 的 可 实施 性 。 
3.2 HY2-SCAT 仿真 

为 了 获取 重 构图 像 ,首先 应 获得 散射 计 每 次 测 
量 的 后 向 散射 系数 > 和 孔径 画 数 刀 和 矩阵 ,在 真实 情 
况 下 ,散射 计 的 测量 值 就 是 后 向 散射 系数 ,但 是 在 仿 
真 时 ,应 该 首先 模拟 散射 计 的 扫描 过 程 ,通过 计算 得 
出 每 一 次 扫描 后 向 散射 系数 > 的 值 ,并 同时 得 出 每 
一 次 扫描 相应 的 刀 和 矩阵。 

仿真 区 域 大 小 设置 为 400 kmX 400 km。 设 重 
构图 像 的 像 元 大 小 为 4, Xd,, 将 地 面 区 域 离散 化 为 
N, XN, 个 网 格 , 并 满足 N. Xd, =L., N, X d, = 
L. ,此 次 重 构 选择 的 区 为 L, 二 L.=400 km, BE 
行 下 点 以 6. 4 km/s 的 速度 自 坐 标 原点 治 y 轴 正 向 
运动 ,正方 形 仿真 区 域 在 坐标 系 第 一 象限 。 天 线 每 
隔 5. 23 ms 发 射 一 个 脉冲 ,脉冲 长 度 为 1.5 ms, 内 


(a) 原始 图 像 


1 
Fig. 1 
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外 波束 间隔 发 射 。 首 先 计算 每 个 脉冲 照射 到 地 面 中 
心 点 的 位 置 和 视 轴 方向 ;然后 ,通过 地 面 像 元 的 坐标 
判断 它 是 否 在 这 次 照射 的 足 印 椭圆 内 ,得 到 相应 的 
h 矩阵 ;最 后 通过 式 (1) 计 算 每 次 脉冲 照射 到 的 足 印 
的 后 向 散射 系数 。 

取 = 0 时 对 应 的 卫星 地 面 坐 标 为 (0,0), 天线 
方位 角 为 0 ,卫星 治 y 轴 正 向 运动 。 对 任意 t 时刻， 
卫星 照射 地 面 的 观测 几何 由 下 式 表 达 : 


t= N puse X 5. 23X107” (11) 
Sat_location=t X vex (12) 
Pan =tX 2/3. 6 (13) 
footcenter_rz=R X cos (pan) d4) 


footcenter_y =R X sin(@m)+Sat_location (15) 
其 中 ;Ni 表示 脉冲 的 个 数 ,Sat_location 表示 上 时 
刻 卫 星 的 位 置 ,vw 表示 卫星 的 对 地 速度 , y, 表示 天 
线 的 旋转 角度 ,R 表示 波束 的 地 面 斜 距 , 外 波束 为 
872.6 km, 内 波束 为 698 km, (footcenter_zx, footcenter 
_J) 表 示 足 印 中 心 的 坐标 。 通 过 上 述 公 式 , 可 以 得 到 
每 一 个 脉冲 照射 到 地 面 的 足 印 的 中 心 坐标 ,通过 足 印 
的 大 小 则 可 以 判断 地 面 被 照射 到 的 区 域 ,从 而 确定 及 
和 矩阵 。 设 在 足 印 内 的 点 刀 和 矩阵 值 为 1; 不 在 足 印 内 的 
点 尹 矩 阵 值 为 0。 通过 计算 像 元 的 中 心 坐 标 到 足 印 


像 元 是 否 在 足 印 内 。 如 果 小 于 或 等 于 长 轴 , 则 在 足 印 
内 ,如 果 大 于 长 轴 , 则 在 足 印 外 ,公式 表达 如 下 : 


( Ts di +d, < 2a 
y “ina +d, > 2a 
(16) 
其 中 :di d: 是 像 元 中 心 到 足 印 椭圆 两 焦点 的 距离 ， 


a 是 足 印 椭圆 的 半 长 轴 长 度 。 
本 文采 用 图 像 处 理 的 经 典 图 像 fruits, 如 图 1(a)， 
作为 仿真 图 像 ,测试 在 不 同 像 元 大 小 下 AART、 


(b) 仿真 测量 图 像 


仿真 图 像 


Original image 


498 E R 


MART 和 SIR 算法 的 重 构 效果 。 
3.3 AART、MART.、SIR 算法 无 噪 仿真 
3.3.1 无 噪 仿真 结果 
从 图 2 一 4 可 以 看 出 ,在 无 噪声 的 情况 下 ， 
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AART, MART, SIR 算法 都 能 很 好 地 重 构图 像 , 重 
构图 比 仿真 的 散射 计 测 量 图 像 细节 更 丰富 ,效果 更 
清晰 ,分 辨 率 更 高 。 随 着 像 元 减 小 , 重 构图 像 的 纹理 
更 加 接近 真实 图 像 。 


(a) 像 元 大 小 2 kmX2 km 


(b) 像 元 大 小 6 km X 6 km 


(c) 像 元 大 小 10 kmX10 km 


2 AART 无 噪 仿真 结果 图 
Fig.2 Simulation results of AART without noise 


(a) 像 元 大 小 2 kmX2 km 


(b) 像 元 大 小 6kmX6 km 


(c) 像 元 大 小 10 kmX10 km 


3 MART 无 噪 仿真 结果 图 
Fig.3 Simulation results of MART without noise 


(a) 像 元 大 小 2 kmX2 km 


(b) 像 元 大 小 6 kmX 6 km 


(© 像 元 大 小 10 km X 10 km 


4 SIR 无 噪 仿真 结果 图 


Fig.4 Simulation results of SIR without noise 


3.3.2 无 噪 情况 下 3 种 算法 的 对 比 
从 表 1 可 知 ,无 噪声 情况 下 AART.MART 和 
SIR 3 种 算法 在 相同 像 元 尺寸 的 条 件 下 ,相关 系数 


和 均 方差 基本 相同 ,都 能 很 好 地 重 构图 像 ,SIR 算法 
MRE =F AART 和 MART 算法 。 


如 图 5 PAM WIS SAM KAMA 
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次 数 的 变化 可 以 得 出 ,在 无 噪声 的 情况 下 AART, 
MART 和 SIR 3 种 算法 在 迭代 若干 步 后 均 可 以 达 
到 收敛 ,AART 算法 像 元 的 更 新 项 趋 于 0, MART 
和 SIR 算法 像 元 的 更 新 项 趋 于 1, 完 成 对 图 像 的 重 
构 。AART 算法 收敛 速度 最 快 ,可 以 以 最 少 的 迭代 
次 数 完成 对 图 像 的 重 构 , 同 时 重 构 后 的 图 像 和 原 图 
像 的 相关 系数 最 高 , 重 构 效果 最 好 ; MART 算法 收 
GR ERF AART 算法 ,但 也 能 较 快 地 完成 收敛 ， 
重 构 后 的 图 像 和 原 图 像 的 相关 性 比 AART 算法 稍 
差 ,也 能 较 好 地 逼近 原 图 像 ;SIR 算法 的 收敛 速度 最 
慢 , 相 关系 数 最 低 ,但 是 在 收 你 之 后 也 能 较 好 地 逼近 
原 图 像 ,和 AART、MART 算法 的 性 能 相差 不 多 。 
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#1 无 噪 情况 下 AART、MART、SIR 算法 性 能 比较 


Table 1 Comparisons of AART,MART,SIR without noise 
算法 名 称 像 元 大 小 相关 系数 均 方 差 
AART 0. 96 11.9 
MART 2 kmX2 km 0.96 11,2 
SIR 0. 96 12. 4 
AART 0. 99 7.0 
MART 6 kmX6 km 0. 98 7.4 
SIR 0. 98 8.3 
AART 0. 99 2.5 
MART 10 kmX10 km 0. 99 ca 
SIR 0.99 5. 0 


3) 20 40 60 80 100 120 140 160 180 200 
ERP 
(a) AREK REE RY LRA 
—— AART 


通过 对 比 各 算法 均 方 差 随 迭 代 次 数 的 变化 可 以 
得 出 ,在 无 噪声 情况 下 AART 、MART SIR 3 种 算 
法 都 可 以 很 好 地 逼近 原 图 像 。AART 算法 收敛 后 
的 均 方差 最 小 ,能 够 最 大 程度 地 接近 原 图 像 ; 
MART 虽然 收敛 较 慢 ,但 在 迭代 150 次 后 , 均 方 差 
5 AART 相差 很 小 , 仍 能 很 好 地 逼近 原 图 ;SIR 算 
法 逼近 速度 最 慢 , 收 剑 后 的 均 方 差 比 AART 和 
MART 略 大 ,但 是 仍 可 以 很 好 地 逼近 原 图 像 ; 在 3 
种 算法 迭代 200 次 后 , 均 方差 的 差别 很 小 ,算法 性 能 
相差 不 大 。 

3.4 HY2-SCAT 噪声 模拟 

散射 计 的 测量 误差 主要 来 自 于 回 波 信号 的 测量 
精度 、 内 定 标 精度 和 天 线 波 束 指向 偏差 。 

散射 计 系统 中 一 般 采用 o 测量 值 的 相对 标准 偏 
= K, 来 评估 回 波 信 号 的 测量 精度 1,K, 决定 于 系 
统 的 信 噪 比 和 进行 平均 的 独立 测量 样本 的 数目 。 

定义 式 如 下 : 


ey 20 40 60 80 100 120 140 160 180 200 
迭代 步骤 
(b) 各 算法 均 方差 随 兴 代 次 数 的 变化 关系 
MART —-—- SIR 
图 5 无 噪 条 件 下 各 算法 性 能 随 迭 代 次 数 变化 关系 
Fig.5 Comparisons of AART,MART,SIR without noise 


Kp= a tons (1 Ea (a) ] a 
HRAN o 测量 结果 可 以 用 式 (19) FAR 
得 co 
二 (19) 
其 中 :7T, 是 接收 时 间 ,B, 是 接收 机 信号 带宽 ,B, 是 
接收 机 噪声 带宽 ,SNR 是 散射 计 信 噪 比 , 信 噪 比 可 
以 由 下 式 计 算得 出 : 


SNR = Et 


P, 
其 中 :P, 表示 系统 接收 功率 ,PP, 表示 系统 噪声 功率 ， 
其 中 P, 可 以 用 式 (21) 计算 得 出 : 
_ PGR AS 
" Cr) R'Lr 
其 中 :P, 表示 系统 发 射 功 率 ,G 表示 天 线 增益 ,) 表 


(20) 


(21) 


500 iz R 


示 雷 达 波 长 ,A 表示 足 印 面积 ,R 表示 斜 距 ,Lr 表示 
系统 差 损 。 

P, = kg e (np—1) * Try © B, (22) 
其 中 :ks 表示 波 尔 兹 曼 常 数 ,nj 表示 接收 机 噪声 系 
数 , 取 值 为 5 dB, Ts 表示 环境 温度 , 取 290 KU, 
3.5 AART、MART.、SIR 有 噪 仿真 
3.5.1 AART 有 噪 仿真 结果 

如 图 6 所 示 ,在 有 噪声 情况 下 AART 算法 的 重 
构 能 力 有 限 ,图 像 的 细节 可 以 恢复 一 部 分 ,但 是 , 因 
为 每 个 像 元 的 更 新 项 是 测量 值 和 前 向 估计 的 差 值 ， 


(a) 像 元 大 小 2 kmX2 km 
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(b) 像 元 大 小 6 km X 6 km 
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在 算法 运行 过 程 中 会 对 图 像 中 的 高 频 分 量 进行 放 
大 ,所 以 图 像 中 存在 很 强 的 点 状 噪声 。 
3.5.2 MART HEAR 

如 图 7 所 示 , 在 有 噪声 的 情况 下 ,MART 算法 
的 性 能 比 AART 算法 好 。 因 为 MART 的 和 迭代 更 新 
项 是 乘 性 的 ,每 一 个 像 元 的 更 新 项 是 测量 值 与 前 向 
估计 的 比值 ,再 与 当前 估计 值 相 乘 ,完成 对 像 元 的 更 
新 。 在 有 噪声 的 条 件 下 , 像 元 的 更 新 项 逐渐 收敛 .在 
单位 值 1 附近 震荡 ,所 以 算法 对 高 频 分 量 的 放大 作 
用 有 限 , 有 一 定 抑制 噪声 的 能 


(c) 像 元 大 小 10 km X 10 km 


6 AART 有 噪 仿真 结果 图 
Fig.6 Simulation results of AART with noise 


(a) 像 元 大 小 2 kmX2 km 


(b) 像 元 大 小 6 km X6 km 


(c) 像 元 大 小 10 km X10 km 


7 MART 有 噪 仿真 结果 图 
Fig.7 Simulation results of MART with noise 


3.5.3 SIR 有 噪 仿真 结果 

如 图 8 所 示 , 在 有 噪声 情况 下 SIR 算法 的 复原 
能 力 最 强 ,噪声 影响 很 小 ,图 片 细节 清晰 , 重 构 效 果 
好 。 这 是 因为 ,每 个 像 元 的 更 新 项 不 再 是 简单 的 测 
量 值 与 前 向 估计 的 比值 ,而 是 通过 对 比例 因子 di 大 
小 的 判断 改变 更 新 项 。 在 存在 噪声 的 情况 下 ,比例 
因子 不 会 最 终 收敛 到 1 ,但 是 因为 算法 根据 比例 因 
子 的 大 小 提供 了 不 同 的 非 线 性 更 新 项 ,这 就 限制 了 
对 前 向 估计 和 测量 值 之 间 的 大 的 偏差 的 更 新 ,从 而 


限制 了 最 终 估 计 值 的 大 小 ,对 噪声 起 到 了 限制 作用 ， 
拥有 很 强 的 噪声 抑制 作用 。 
3.5.4 ARIEN TF 3 种 算法 的 对 比 

MR 2 可 知 ,在 有 噪声 的 情况 下 AART 算法 的 
相关 系数 和 均 方差 最 差 ,MART 算法 的 相关 系数 和 
均 方差 优 于 AART 算法 ,SIR 算法 的 相关 系数 和 均 
方差 最 好 。 这 说 明 SIR 算法 对 噪声 的 抑制 效果 最 
好 ,MART 有 一 定 的 抑制 噪声 效果 ,而 AART 对 噪 
声 几 乎 没有 抑制 效果 。 


paN 
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(a) 像 元 大 小 2 kmX2 km 


Fig. 8 
表 2 有 了 噪 情况 下 AART、MART、SIR 算法 性 能 比较 
Table 2 Comparisons of AART, MART,SIR with noise 
算法 名 称 像 元 大 小 相关 系数 均 方 差 
AART 0. 89 218 
MART 2 kmX2 km 0.93 6.8 
SIR 0.95 3.6 
AART 0. 92 8.6 
MART 6 kmX6 km 0:95 3 
SIR 0.97 0.0 
AART 0. 95 4.2 
MART 10 kmX10 km 0.97 0.9 
SIR 0.99 Lace 


如 图 9 所 示 WIS SAM KAMAN 
次 数 的 变化 可 以 得 出 ,在 有 噪声 的 情况 下 AART, 
MART 和 SIR 算法 随 着 迭代 步骤 的 增加 ,最 后 不 能 


收敛 , 即 AART 算法 的 更 新 项 无 法 收敛 至 0, 而 
MART, SIR 算法 的 更 新 项 无 法 收敛 至 1, 始终 存在 


(b) 像 元 大 小 6 kmX6km 
8 SIR Blk 


(c) 像 元 大 小 10 km X10 km 


Bet 


aso 


RE 


Simulation results of SIR with noise 


=m. AART RAHI H X A A E E A NR A Iie 
加 , 先 迅 速 上 升 随后 快速 下 降 ; 相 应 地 均 方 差 随 着 和 迭 
代 次 数 的 增加 , 先 迅速 下 降 随 后 快速 上 升 。MART 
算法 的 相关 系数 随 迭 代 次 数 的 增加 先 迅速 上 升 后 快 
速 下 降 , 但 是 上 升 速率 和 下 降 速 度 和 幅度 都 低 于 
AART 算法 ;相应 地 , 均 方差 随 着 迭代 次 数 的 增加 
先 下 降 后 上 升 , 上 升 和 下 降 的 速度 和 幅度 都 小 于 
AART 算法 。SIR 算法 的 相关 系数 随 着 迭代 次 数 
的 增加 上 升 , 上 升 速度 低 于 AART、MART 算法 ， 
但 是 在 上 升 至 峰值 后 , 随 着 和 迭代 次 数 的 增加 ,下 降幅 
度 很 小 ,基本 保持 在 很 高 的 水 平 ;相应 地 , 均 方 差 随 
着 和 迭代 次 数 的 增加 , 先 快速 下 降 , 随 后 有 小 幅度 上 
升 ,但 是 基本 保持 在 一 个 很 低 的 水 平 。 造成 
AART、MART 和 SIR 算法 在 有 了 噪 情 况 下 ,相关 系 
数 和 均 方差 变化 的 原因 在 于 :图 像 的 误差 由 两 部 分 
构成 ,一 部 分 是 重 构 误 差 , 另 一 部 分 是 噪声 放大 。 算 
法 会 降低 重 构 误差 ,但 是 会 将 存在 于 图 片 中 的 噪声 


0 100 200 300 400 
(a) 相关 系数 
— AART 


500 600 700 800 900 1000 


105700 200 300 400 
(b) AAE KEL 
=== SIR 
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9 有 噪 条 件 下 各 算法 相关 系数 和 均 方 差 随和 迭代 次 数 变化 关系 
Fig.9 Comparisons of AART,MART,SIR with noise 
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声 的 放大 ,所 以 相关 系数 上 升 、 均 方 误差 下 降 , 逐渐 
逼近 算法 性 能 的 最 佳 值 ;但 随 着 迭代 次 数 的 增加 , 算 
法 对 于 重 构 误差 的 降低 小 于 对 噪声 的 放大 ,所 以 出 
现 了 相关 系数 下 降 和 均 方 误差 增加 。 从 图 9 可 以 看 
出 ,AART 算法 对 噪声 的 抑制 最 差 , 所 以 在 最 佳 值 
过 后 相关 系数 和 均 方差 严重 恶化 ;MART 算法 对 噪 
声 有 一 定 的 抑制 能 力 , 在 最 佳 值 过 后 相关 系数 和 均 
方差 的 恶化 程度 小 于 AART 算法 ;SIR 算法 对 噪声 
的 抑制 最 好 ,在 最 佳 值 过 后 ,相关 系数 只 是 小 范围 下 
降 , 均 方差 只 有 小 范围 的 增加 , 总体 处 在 很 好 的 水 
平 。 所 以 ,在 有 噪声 的 情况 下 SIR 算法 的 性 能 最 
好 。 同 时 ,在 有 噪声 的 情况 下 AART, MART 和 


120 121 122 123 124 125 126 
经 度 /* 


(a) MART 重 构图 
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SIR 3 种 算法 都 不 能 最 终 收 敛 ,图 像 重 构 的 效果 和 与 
迭代 次 数 相关 ,存在 迭代 次 数 的 最 佳 值 , 所 以 在 以 后 
进行 真实 数据 处 理 时 ,要 特别 注意 迄 代 次 数 ,使 算法 
性 能 发 挥 到 峰值 , 重 构图 像 效 果 更 好 。 


4 HY2-SCAT 散射 测量 数据 图 像 重 构 


4.1 HY2-SCAT 数据 
本 文采 用 HY2-SCAT 散射 计 LIB 级 数据 , 提 
取 位 置 区 域 的 经 度 、 纬 度 .方位 角 、 天 线 方位 角 和 后 
向 散射 系数 信息 ,通过 3 种 算法 重 构 目 标 区 域 图 像 。 
本 文选 取 120 °~126 °E, 10 °~ 16 °N 的 区 域 作 为 
重 构 目 标 区 域 ,目标 区 域 为 菲律宾 部 分 海岛 ,散射 计 
原始 后 向 散射 系数 所 构成 的 散 点 图 如 图 10(a) 所 示 。 


(b) AART 重 构图 


(b) SIR 重 构图 


图 10 HY2-SCAT 真 实数 据 重 构 结 


Fig. 10 Image reconstruction of HY2-SCAT data 


从 图 中 可 知 ,散射 计 后 向 散射 系数 直接 做 出 的 图 无 
法 细致 区 分 出 陆地 ,海洋 和 海岸 线 ,不 能 满足 对 陆地 
的 观测 需求 。 
4.2 HY2-SCAT 真实 数据 重 构 结果 

从 图 10 可 知 ,AART、MART 和 SIR 3 种 算法 


都 可 以 从 散射 计 的 后 向 散射 系数 测量 值 中 较 好 地 恢 
复出 地 貌 特征 ,可 以 较为 细致 地 区 分 出 陆地 海洋 和 
海岸 线 。 但 是 3 种 算法 对 噪声 的 抑制 作用 有 着 明显 
区 别 。AART 算法 对 噪声 有 很 明显 的 放大 作用 ,从 
图 中 可 以 看 到 在 陆地 部 分 有 很 多 的 高 频 分 量 ,影响 
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对 陆地 的 观测 ;MART 算法 对 噪声 有 一 定 的 抑制 作 
用 ,从 图 中 可 以 看 到 陆地 上 的 高 频 分 量 明 显 少 于 
AART 算法 ,但 仍 有 较 多 的 点 状 噪声 ; SIR 算法 对 
噪声 的 抑制 作用 最 好 ,从 图 中 基本 无 法 看 到 被 放大 
的 噪声 高 频 分 量 , 图 像 细 节 丰 富 、 轮 廓 清晰 ,噪声 抑 
制 效 果 好 。 


5 结 语 


传统 的 星 载 微 波 散 射 计 分 辨 率 较 低 ,多 用 于 测 
量 海洋 风 场 。 随 着 地 球 探测 需求 的 不 断 增 加 , 星 载 
散射 计 开始 进行 陆地 ,海洋 、 冰 盖 和 植被 的 测量 , 传 
统 星 载 散射 计 分 辨 率 精度 不 能 满足 要 求 , 所 以 需 
提高 星 载 散射 计 的 分 辨 率 。 

本 文采 用 图 像 重 构 算法 对 散射 计 测 量 所 得 后 向 
散射 系数 进行 高 分 辨 处 理 ,仿真 了 AART、MART 
Al SIR 3 种 算法 在 无 噪声 和 有 了 噪声 的 情况 下 图 像 重 
构 ,并 得 出 在 无 噪声 情况 下 ,3 种 算法 都 可 以 很 好 地 
进行 图 像 重 构 , 重 构 出 的 图 像 相 比 仿真 的 测量 图 像 
轮廓 更 清晰 细节 更 丰富 的 结论 。 在 有 噪声 的 情况 
下 ,SIR 算法 拥有 最 强 的 噪声 抑制 能 力 , 重 构 出 的 图 
像 边 缘 清晰 、 细 节 丰 富 ,而 AART、MART 算法 对 
噪声 的 抑制 效果 较 差 ,有 非常 明显 的 点 状 噪声 。 

在 仿真 的 基础 上 进行 了 HY2-SCAT 测量 数据 
的 重 构 ,得 出 了 与 仿真 结果 一 致 的 结论 。AART、 
MART 和 SIR 3 种 算法 都 可 以 较 好 地 重 构 图 像 ,但 
是 SIR 算法 有 较 强 的 噪声 抑制 作用 , 重 构 的 图 像 基 
本 不 存在 点 状 噪 声 , 对 陆地 的 刻画 更 清晰 。 重 构 后 
的 图 像 可 以 应 用 于 植被 . 冰 盖 和 台风 观测 等 地 球 地 
理学 研究 领域 。 


致谢 :感谢 国家 卫星 海洋 应 用 中 心 提供 的 海洋 二 号 
卫星 微波 散射 计 LIB 数据 ! 
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DEM were extracted based on images captured by UAV platform by means of splicing stereo pairs, and 
were compared with GPS measured data, ASTER GDEM data and SRTM data. The result shows that, the 
DEM extracted based on UAV images has the closest difference with GPS measured result (RMSE = 
8. 96) ,which has higher precision than ASTER GDEM (RMSE=13. 68) and SRTM (RMSE=11. 81) da- 
ta; The closest result( RMES= 1. 813) of max tree heights and max elevation differences in each sample 
plot, shows that UAV DEM has an ability to reflect more hierarchical information between canopy and 
ground, which shows a greater potential in the construction of DEM in mountain plantation landscape. 
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Abstract; High-resolution reconstruction techniques are investigated for Hai Yang-2 Scatterometer (HY2- 
SCAT) backscatter measurements for a variety of application purposes including land and polar ice detec- 
tion and the studies of cryophere and biosphere. Three resolution enhancement algorithms, including Addi 
tive Algebraic Reconstruction Technique (AART), Multiplicative Algebraic Reconstruction Technique 
(MART) and Scatterometer Image Reconstruction (SIR),are compared by simulation approach. The re- 
sults that SIR has the best performance in noise suppression and resolution enhancement of the backscatter 
measurements than AART and MART algorithms. This study is verified using HY2-SCAT LeveF1B 
(LIB) backscatter measurements. The o’ measurements over some sites with particular backscattering 
characteristics(e. g. islands) are reconstructed. 
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